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In addition to the generation of waves, a wind produces a mass transport in a body of 
water resulting in the lowering of the level at the windward side and rise at the leeward 
side, which is called wind tide or set-up. Two effects of the wind are involved: the surface 
traction on the water, and the form resistance of the waves. This paper presents the 
theoretical background of the subject and experimental results. It was discovered that 
formation of waves in an experimental channel could be inhibited by adding small amounts 
of soap or detergent to the water. This made it possible to study the surface traction 
effect separately. The effect was studied for both laminar and turbulent motion of the drift 
and gravity currents produced by the wind. The set-up computed from measured wind and 
water surface velocities agreed with theory. Regardless of the flow regime, the set-up 
was unaffected by viscosity and independent of the depth-length ratio of the channel. 
The additional set-up due to presence of waves could be correlated only by introducing a 
characteristic velocity. This additional set-up appears to vary as the square root of the 
depth-length ratio. The relation of the characteristic velocity to the critical velocity for 
wind generation is discussed. The derived empirical formula for set-up is compared with 
observations in Lake Erie. 



I. Introduction 

Early in the course of a comprehensive investiga- 
tion of density currents now under way, the question 
arose as to what phenomena might have a probable 
bearing, no matter how remote, on the general prob- 
lem of density currents and especially on the stability 
of the interface between two fluids of different den- 
sities. Among such phenomena are the disturbances 
that may be introduced by wind-generated surface 
waves on the interface between two layers of dif- 
ferent densities such as occur in lakes and in the 
lower reaches of rivers emptying into the ocean. 

At least two ways in which surface wind waves 
may affect the interface may be imagined. The first 
possibility is that internal waves will be generated 
as a response to the surface waves. This, however, 
is very unlikely to happen. It is well known that a 
two-layer system of liquids with a free surface is a 
system of two degrees of freedom. To one degree 
of freedom there corresponds the ordinary waves of 
a homogeneous liquid, and to the other the internal 
waves of the interface. The velocities of propaga- 
tion of these two types of waves are unequal, and 
the particle motions are significantly different from 
each other. In ordinary waves of a homogeneous 
liquid, the particle motions along a given vertical 
line are in the same direction from the bottom to the 
free surface. With internal waves on the other 
hand, the particle motions are oppositely directed in 
the two layers of different densities. It follows that 
whatever may be the agency that causes the ordinary 
type of waves in a homogeneous liquid, the same 
agency cannot produce internal waves if the homog- 
eneous system is replaced by a system of two layers 
of liquids. Wind waves are surface waves of ordi- 
nary type and therefore will not excite internal 
waves. If internal waves are to be created, they 
must be produced by displacements of the internal 



surface that will leave the free surface practically 
smooth. 

The second possible effect of surface wind waves 
arises from their tendency to produce turbulence in 
the region occupied by the interface, which causes 
mixing of the two liquids. It is conceivable that 
this mixing is brought about in two ways. First, 
wind-generated waves are associated inherently with 
turbulence, and a progressive downward movement 
of turbulence may reach the level of the interface. 
If the density difference at the two sides of the inter- 
face is small, and if the turbulence fluctuations are 
sufficiently strong, a gradual mixing of the liquids 
obviously will take place. 

The second and more significant way in which 
mixing at the interface might occur may be illus- 
trated by the case of a lake subjected to wind wave 
where there is a considerable wind tide or set-up 
present. If there are two layers of liquids of different 
densities in the lake, the set-up will produce circula- 
tions of opposite directions in the two layers and 
thereby will cause mixing. If the density differences 
are small, a similar effect may also be expected to 
occur in the case of salt intrusion in a river when the 
wind is blowing in the upstream direction. The 
smaller the river velocity, the larger will be the mixing 
due to wind effects. 

Before considering the effect on interfacial sta- 
bility of surface disturbances caused by wind, it is 
necessary that sufficiently full information be ob- 
tained on the dynamical and geometrical character- 
istics of the wind waves and especially of wind waves 
in shallow and long channels. There appear to be 
no complete investigations, either theoretical or 
experimental, regarding such waves. Having in 
mind the excellent theoretical work of Sverdrup and 
Munk [1] 1 on wind waves and swells originating in 
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deep waters, it was thought that an experimental 
investigation relating to shallow waters would be a 
relatively simple matter. It was soon evident that 
there are many difficulties in the study of wind waves 
in a laboratory channel with shallow water that are 
absent in the case relating to deep water. An ex- 
tensive program of investigation is practically un- 
avoidable. As the subject is of much wider interest 
and practical importance than merely in its relation 
to density currents, a separate project was initiated. 
For convenience this project has been divided into 
three parts: (1) wind tides or set-up in small channels, 
(2) the dissipation characteristics of oscillatory pro- 
gressive waves in long shallow channels, and (3) the 
growth of wind waves in height and length in long 
shallow channels. 

The present paper deals with the first subject, 
that is, wind tides or set-up in small closed channels. 
This subject is of considerable practical importance 
in relation to wind effects on long shallow lakes and 
the feasibility of model tests. The problem of Lake 
Okeechobee falls into this category. Although a 
final solution of the question of model testing and 
of valid transference equations is not reached, never- 
theless the results obtained represent a method of 
approach toward this solution. 

II. Theoretical Developments 
1. Generalized Surface Equation 

The generalized surface equation is the differential 
equation involving the surface displacement of a 
limited body of water, the free surface being sub- 
jected to a tangential stress. We shall consider a 
channel closed at its ends with a constant width 
throughout and the bottom horizontal. The surface 
equation can be established very conveniently by 
considering the general equations of motion of the 
Navier type. 

Selecting the origin of the axes at the bottom, 
with x along the channel axis and z vertical and 
drawn upwards as in figure 1, the equations of 
motion are: 
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Figure 1. Coordinate axes and stresses. 



Here u and w are the velocity components parallel 
to the x and z axes, respectively, p is the pressure, 
g the acceleration of gravity, and p the density of the 
liquid. The quantities N u N3, T 2 , are the tractional 
stresses as shown in figure 1. In the notation of 
Lamb, N 1 =p+p xx , N s =p+p zz , and T 2 =p X2 . 

The equations in the form given are too general 
and too difficult to integrate. In the case of a long 
shallow channel some simplifying assumptions may 
be adopted. These assumptions are commonly re- 
sorted to in solutions of boundary layer flow prob- 
lems and are well known. As the surface curvature 
is small, the vertical acceleration and the normal 
traction N s will be neglected. As the surface dis- 
placement is small, the normal traction iVi will be 
neglected. The variation of T 2 with respect to x 
is negligible since the channel is very long. Ac- 
cordingly, the simplified forms of the equations of 
motion, writing r for T 2 are 
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It is to be noted that one consequence of the assump- 
tions is to regard the pressures as hydrostatic. For 
the long central portion of the channel, these forms 
of the equations are sufficiently valid. They fail, 
however, to represent the conditions of flow at the 
ends of the channel. 

Denoting the undisturbed depth of the water by 
H 1 as in figure 2, the surface displacement as meas- 
ured from the initial level by h, and the atmospheric 
pressure by p a , eq 4 yields 



p=pg(H+h — z)+p a . 

Substituting in eq 3, 
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Now, since the liquid is incompressible, 
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The rule 



bu b f s bw 
w ^-=^— (wu)—u -5— j 
bz bz bz 



combined with eq 7 yields 



bu d , v . bu 



bz bz 

Substituting this in eq 6, 



du du 2 _d^ , w _ bh 1 bp a 1 br 



(8) 



(9) 



(10) 



Multiplying the two sides of the equation by dz, 
integrating with respect to z between the limits and 
(H+h), and noting that w=0 at the bottom where 

2=0, 
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where w s and ^ s are the velocity components at the 
free surface, r s and — t are the tangential stresses 
at the free surface and at the bottom (see fig. 2). 
We may express w s in terms of v 8 by considering the 
surface kinematic condition that a particle at the 
free surface remains on the surface. Symbolically, 
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Introducing this in eq 11 and remembering Leib- 
nitz's rule for differentiation under the integral sign, 
we have 
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u 2 dz = 






(13) 



This is the generalized differential equation of the 
surface for a long body of water subjected to the 
surface tractions t s and — r . The terms on the 
left-hand side are the acceleration terms. 

If the conditions are steady, as when the wind has 
acted with uniform intensity over the body of water 
for a long time, the equation reduces to 
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Under certain assumptions, to be discussed later, 
the acceleration term on the left-hand side of the 
equation may be neglected. Thus, the surface 
equation reduces to 
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From this point on it will be supposed either that 
the atmospheric pressure or, more properly stated, 
the pressure on the free surface, does not vary with 
x; or that if it does vary, a correction is made so that 
h will denote the surface elevation after due account 
is made of the pressure effect. With this under- 
standing, the basic simplified surface equation is 
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A further simplification, based on assumptions 
relating to the relative values of r and r s , is needed. 
The flow in the channel is a combustion of a drift 
current at the surface and a gravity current at the 
bottom. The drift current is induced by the wind. 
The gravity current is induced by the inclination 
of the upper free surface, its direction being opposite 
to that of the drift current. 

Denoting the mean velocity of the gravity current 
by u , the resistance of the channel bottom may be 
written as 



Tq~- 
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Similarly, if V is the wind velocity and p a the density 
of the air, the surface traction may be expressed as 



s = k sPa V 2 . 



Considering the ratio 
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it is not obvious that r is negligible with respect to 
r s . The velocity u may be small with respect to V, 
but since p is about one thousand times larger than 
p a , the product pu% may be a measurable portion 
of p a V 2 . Besides one must also consider the relative 
values of k s and k Qj and at the moment it cannot 
be said that k is smaller than k s . For these reasons, 
it is better not to neglect r with respec to r„ but 
to incorporate it with r 8 . Writing 



TO 
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eq 17 becomes 
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If a computation is made on the basis of pure 
viscous laminar motion, and this will be done in one 
of the following sections, it will be seen that n=1.5. 
The determination of n for turbulent motion requires 
some form of turbulent theory for the case under 
consideration. The value of n will then depend on 
the particular theory adopted. Hells trom [2], after 
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adopting the Boussinesq theory for turbulent flow, 
finds that n varies between the values 1.15 and 1.80 
for channels of moderate and very large depths. 
Temporarily we shall take n=1.25 for turbulent 
flow. 

2. Surface Equation for Laminar Flow 

From the viewpoint of theory it is important to 
consider the case of laminar flow. The chief interest, 
of course, rests on the fact that for this case a solu- 
tion is easily available. In laboratory experiments, 
as will be shown later, this solution is of great value 
when the depths are small. 

In viscous flow of the type envisaged 






(23) 



Substituting in eq 10, neglecting inertia terms, sup- 
posing the flow to be steady, the pressure over the 
surface uniform, and the surface displacements small 
with respect to the depth of water, there results 
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If the velocities are expressed in terms of the surface 
velocity, u s , and the distances from the bottom in 
terms of the initial liquid depth, H, this equation may 
be made dimensionless. Thus, writing 
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eq 24 becomes 
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which is'subject to the boundary conditions 
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The first condition states that the velocity at the 
bottom is zero, the second that velocities are meas- 
ured in terms of the surface velocity, and the third 
that the integrated flow in a vertical cross section 
vanishes. 



The desired solution subject to the boundary con- 
ditions is 

0=3f 2 -2f, (29) 

and 

a=6. (30) 



Equation 29 gives the law of velocity distribution. 
In the horizontal plane at a distance of II/3 below 
the liquid surface the velocities vanish. All the 
particles above this plane move in the direction 
of the wind, and the particles below this plane 
move in the opposite direction. The former repre- 
sents the drift current, and the latter the gravity 
current. Equation 30 gives the surface equation; 
for, using eq 26, 
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which relates the surface inclination to the depth 
of liquid, the surface velocity, and the kinematic 
viscosity of the liquid in the channel. 
Now, from eq 23 and 25, 
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The surface shear and the shear at the bottom 
correspond to f = l and f = 0, respectively, and we 
have for these shears 
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(34) 



that is, the surface shear is twice as large as the shear 
at bottom. The surface equation may be written 
also as 



— =1 5 T * 
dx pgH 



(35) 



It is apparent that ignoring the bottom shear results 
in a 50-percent error. 

III. Characteristics of the Experimental 
Channel 

1. Apparatus 

The apparatus, shown diagrammatically in figure 
3, consisted of three parts: the entrance assemblage, 
the experimental channel, and the terminal assem- 
blage. 

The upstream blower in the entrance assemblage 
was driven by a constant speed motor. The air 
volume, and hence the air velocity in the experimental 
channel, was controlled by a damper at the intake 
end of the blower by means of which the intake area 
could be varied. The control of air velocities by 
this method was very satisfactory. 

At the discharge end of the blower was a rectan- 
gular duct of the same dimensions as the blower 
discharge opening connected to a duct converging 
in plan to the width of the experimental channel. 
This was followed by a Venturi nozzle of constant 
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Figure 3. Diagram of experimental apparatus. 
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width and converging in the vertical dimension to 
half the height of the experimental channel. Pie- 
zometer openings in the side walls of the Venturi 
nozzle were connected to an inclined manometer. 
The Venturi was connected by a short duct to the 
upper half of the experimental channel. It was 
thought that the length of the assemblage and the 
lateral and vertical contractions would even out 
most of the irregularities in flow produced by the 
blower. 

The experimental channel, shown in figure 3, was 
a closed conduit of uniform rectangular cross section, 
28.5 cm deep, 11.3 cm wide, and about 20 m long. 
It consisted of 12 sections constructed of wood with 
sides of glass. The glass sections were connected 
by short sections made of wood. Piezometers, 
opening to the inside air at the top and to water 
at the bottom of the channel were installed in the 
walls of the wooden sections. The piezometer 
openings were connected to glass manometer wells. 
The average water levels in the channel were meas- 
ured by the water levels in the manometers. 

It was desired to make tests with different depths 
of water. As it was not practicable to change the 
elevation of the Venturi nozzle, a sloping floor was 
set in the upstream end of the channel. The rate 
of variation of the air cross section was quite small, 
thus assuring a nearly uniform air velocity distri- 
bution vertically. 



In the top of each channel section at the center 
was placed a piezometer opening, and the drop in air 
pressure between any pair of these could be meas- 
ured by connecting them to an inclined manometer. 
During observations of set-up, those two piezometers 
were used that were nearest to the ends of the fetch 
over which measurements were being made. 

The terminal assemblage, also shown in figure 3, 
began at the end of the experimental channel with a 
section expanding in width. A succession of screens 
was placed at the lower end for damping of the 
waves. The screens were spaced about 2 cm apart 
initially and the spacing decreased to about X cm 
at the downstream end. The arrangement proved 
to be quite satisfactory. The damping was gradual 
and effective, there being no appreciable reflection 
of the waves. Downstream from this section was 
a duct sloping slightly upward and expanding to 
the size of the intake opening of the downstream 
blower. This blower was used in addition to the 
upstream one whenever high air velocities were 
required. 

2. Calibration of the Venturi Nozzle 

The mean velocity of the air stream over the 
water in the channel was chosen as the significant 
velocity to be associated with the growth of waves 
or with the development of wind set-up. The depth 
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of the channel being small, the height of the air 
passage over the water decreased continually in the 
leeward direction because of the set-up and of the 
increase in wave heights. Therefor it was not very 
practical to measure the local mean air velocity at 
any chosen cross section. Instead, the air velocity 
was deduced from the velocity at the discharge end 
of the Venturi nozzle. Thus if T 7 i is the velocity 
and A\ the area of the nozzle exit, and T 7 is the mean 
air velocity in the channel for an air passage of 
cross-sectional area A, then V=ViAi/A. The area, 
A , represents the mean value of the area averaged 
along the entire length of the channel. Obviously, 
A is also the area of the air passage cross section 
when the water in the channel is not disturbed. In 
the measurements of set-up, however, the decrease 
of the cross section of the air passage with distance 
was considerable, and hence corrections were devised, 
as will be explained later, to obtain the value of the 
effective wind velocities. 

The differential head reading of the inclined 
manometer was an indication of the velocity at the 
discharge end of the Venturi nozzle. The calibra- 
tion was made in the following manner. The ter- 
minal assemblage at the lower end of the experimental 
channel was removed, and the velocities at the end 
section of the channel were measured. These were 
associated with the simultaneous readings of the 
Venturi manometer. The velocities at the channel 
end were measured by two instruments. One was 
a propeller-type anemometer of current manufacture 
that was calibrated in a wind tunnel by the Aero- 
dynamics Section of the National Bureau of Stand- 
ards. The second instrument was a small Ott cur- 
rent meter. Although this instrument is used 
ordinarily for flowing water and the calibration is 
likewise made in water, it can, nevertheless, be used 
for flowing air provided the calibration of the instru- 
ment is expressed in the form 



N/U=f(Uifi), 



(36) 



where N is the number of revolutions per second 
corresponding to the velocity, U. [3] To be sure 
that this relationship applied to the Ott meter to 
be used, it was calibrated both in water and in air. 
In the range where the calibration results with air 
and with water overlapped, the results fell on the 
same curve. 

The Ott meter, being the smaller of the two, was 
set in six different positions in the channel cross 
section. The other instrument was set in three 
positions on the vertical center line. The average 
of the mean velocities determined in these two ways 
was taken to be the mean velocity of the air stream 
at the end cross section of the channel. Comparing 
the cross-sectional areas, the velocity at the nozzle 
was obtained in relation to the manometer reading. 
Following this, a second calibration was obtained 
after attaching to the end of the main channel an 
additional closed rectangular conduit of a larger 
cross section. In the second channel both instru- 
ments were set in six different positions. 



It was gratifying to note that the Venturi nozzle 
velocities as observed during the calibrations could 
be very well represented by the theoretical expression. 

2 / ><M = Pa y?[l-(^) 2 ], (37) 

where Ah is the reading of the inclined manometer in 
feet of water, and A 2 is the area at the entrance of the 
nozzle. The final calibration chart was made by 
using the above formula and taking pjp = 1.1 6X 10~ 3 , 
a value corresponding to an atmospheric pressure of 
740 mm Hg and an air temperature of 23° C. 



3. Air Velocity Distribution in Channel 

In studies of the growth of waves with distance 
and of the development of wind set-up under natural 
conditions the proper selection of the elevations for 
the measurement of wind velocities is a difficult 
matter. If one could assign an operative roughness 
length, k, independent of the dimensions of waves 
excited by the wind, it would be sufficient to measure 
the wind velocity at a constant elevation, z aj inde- 
pendent of the wind velocity. The investigations of 
Roll [4] concerning the wind field over the seas show T 
that an operative constant roughness in the above 
sense does not exist and the roughness is dependent 
on the wind velocity. This naturally introduces a 
complication if results of investigations in the open 
seas are to be compared with results obtained under 
laboratory conditions. 

Because of the limited channel dimensions there 
is little latitude in selecting the location for the wind 
velocity measurements in laboratory tests. The 
velocity of the wind may be taken either as the 
velocity at the center of the area above the water or 
as the mean velocity all over the area. In the 
present investigation the second choice was preferred. 
This procedure would be free of any objection pro- 
vided the air velocity distributions in the channel 
remain always afflne to each other for varying 
velocities and varying water depths. This matter 
was examined for a water surface free of waves. 

One additional point that needed examination was 
the assumed effectiveness of the expansion caused 
by the sloping bottom at the entrance of the main 
channel (see fig. 3) in establishing uniform air ve- 
locities at the leeward end of the expansion. For 
this purpose two sections were chosen, one at the 
windward end and the other at the leeward end of 
the channel, and vertical air velocity traverses 
were made along the vertical median plane. One 
group of such determinations is shown in figure 4, 
where V represents the velocity at a distance, z a , 
from the water surface, V m the average velocity 
along the vertical, and II a the height of the air 
gap. In this case the water depth was 4.0 cm, the 
air space height 24.5 cm, and the average velocity 
along the vertical 1,008 cm/sec. The velocity 
distributions at the two ends of the experimental 
channel are practically alike. For the same air 
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space similar determinations were carried out for 
two other velocities, namely 700 and 480 cm/sec, 
with results very similar to those shown in figure 4. 
The measurements were repeated for other air gaps, 
namely H a =20.5, 17.5, 14.0 cm, with three air 
velocities not greatly different from those mentioned 
above. It was noted in every case, regardless of the 
air space or the magnitude of the mean air velocity, 
that the median vertical velocity distributions were 
alike in the extreme end cross sections. Thus for 
the comparison of the velocity distributions in each 
air space it should suffice to consider the averages 
from the measurements with a chosen air space. 
Such mean values are shown in figure 5. From the 
curve drawn in the figure we make two deductions. 
First, central wind velocities are about 1.1 times as 
large as the mean air velocity. Second, there are 
no large differences between the air velocities near 
the rigid top of the channel and those near the water 
surface. If anything, the air velocities near the 
water surface are slightly higher, as might be 
expected. 

Next, the horizontal velocity distribution close to 
the water surface was examined. The results of a 
few traverses are shown in figure 6. The plotting is 
V/V m&x versus y/B, where V m&x is the velocity at the 
median plane, V the velocity at a distance, y, from 
one of the vertical walls, and B the width of the 
channel. The liquid in the channel was of large 
viscosity so that for large air velocities the surface 
would be free of waves. Three series of traverses 
were obtained for the surface of the liquid thus un- 
modified by waves corresponding to settings of the 
pitot tube at successive elevations, z a , of 0.4, 0.9, and 
1.7 cm above the water surface. The results are 
shown in the figure by the open symbols. Only one 




Figure 4. Vertical distributions of air velocities at the ends of 
the channel. 

Circles, velocities at windward end; squares, velocities at leeward end. 



measurement of velocities was made with the liquid 
surface covered with waves. In these measurements 
the pitot tube aperture was almost at the level of the 
crests of the waves. The data for these tests are 
shown by black circles. It is important to note that 
for the lateral traverses the distributions of the 
velocities are more uniform than in the vertical 
traverses. The presence of the waves hardly modi- 
fies this comparison. 

IV. Wind Set-Up in the Absence of Waves 
1. Effect of Soapy Waters on the Start of Waves 

It was discovered accidentally during the conduct 
of the routine experiments that the addition of soap 
to the water inhibited the formation of waves by 
wind action. Ordinarily when the water in the 
channel was free of foreign matter, a wind with 
velocity of 300 to 400 cm/sec caused waves to be 
formed. An increase in the wind velocity above this 
value increased the height and length of the waves. 
On the other hand, when soap was introduced in the 
water and the wind velocity was increased gradually 
from zero value, it was found that no matter how 
great the wind velocity, the surface of the liquid re- 
sisted the formation of waves. Even at a wind 
velocity of 1,200 cm/sec, which was the limiting 
capacity of the apparatus, the soapy liquid surface was 
free of waves. This peculiar phenomenon could be 
repeated for all depths of liquid from 4 to 14.5 cm. 
The same condition could also be reproduced by 
introducing a detergent of the type sold for cleansing 
kitchen utensils, instead of soap. 
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Figure 5. Vertical distributions of air velocities at center-line 
of channel. 
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Figure 6. Horizontal distributions of air velocities near liquid 
surface. 



It must not be inferred from the above that in con- 
taminated waters of the kind described all wave 
motions are impossible. If, by a special manipula- 
tion, waves of substantial dimensions are once pro- 
duced, these continue to grow in size and multiply 
in number after traveling long distances under the 
action of the wind. We have had a very significant 
experience demonstrating this fact which deserves a 
brief mention here. 

The soapy water in the channel was about 10 or 
11 cm deep, the two blowers were in action and, owing 
to the large air velocity, the water surface was con- 
siderably inclined, the greater depth, of course, 
being at the leeward end of the channel. The surface 
was smooth and free of waves. The windward edge 
of the water was located on the solid sloping surface 
at the channel entrance. To see what effect a sudden 
reduction of the air velocity would produce, one of 
the blowers was stopped quickly. Immediately a 
surge from the leeward end of the channel started 
to travel toward the windward end. At the instant 
that the surge reached the windward end and the 
water commenced to pile up on the sloping surface 
there, the idle blower was started again. Under the 
impact of the increasing wind two joining intumes- 
censes were seen to develop at the windward water 
edge which moved in the leeward direction. In the 
movement the number of the waves continually 
increased and the crests became higher. Thus a 
limited line of waves, so to speak, was moving down- 
stream, the length of the line increasing continually. 
Waves in the front of the line would be pressing for- 
ward, while at the rear new waves appeared. Finally, 
a long line of waves was filling the last quarter of the 
channel whereas the water in the rest of the channel 
was smooth and undisturbed. This line of waves 
would be seen to enter the wave dampers, and after 



a while the last wave of the line would be lost. 
Following this a new surge started moving toward the 
windward end. The cause of the surge was the sud- 
den decrease of the wave friction concommitant 
with the disappearance of the waves. The surge 
reaching the windward end of the channel, a new 
pair of waves was seen to be formed and the entire 
process of wave group motion described above was 
repeated. The phenomenon appeared to reproduce 
itself periodically and uniformly. Obviously, here 
was a special case involving the principle of group 
waves and their generation under the action of 
wind forces. It is regrettable that, due to the pres- 
sure of the main problems of the investigation, addi- 
tional time could not be devoted to the study of this 
unusual phenomenon. 

Thus the inhibitive action of the soapy water or 
water containing detergent is to be associated with 
water surfaces initially undisturbed and therefore 
free of local elevations occurring over short lengths. 
Under these conditions, no matter how strong the 
wind may be, the action of the flowing air fails to 
start those tremors of the water surface that are 
always present when the water is pure and the for- 
mation of waves is impending. The action of the 
contaminated water cannot be explained. It would 
appear that neither viscosity nor surface tension have 
a bearing on the subject. 

Once the fact became known that it was possible 
to operate with a smooth surface in the channel 
within the full capacity of the blowers even though 
the wind velocities were large, by merely adding 
soap or detergent to the water, it was decided to 
make full use of it. Numerous experiments were 
carried out under these conditions that are considered 
in the subsequent sections. It was felt at the time 
and also verified later that through these experiments 
a better understanding of the main problems of the 
investigation would result. 



2. Formulation of Wind Set-Up for Smooth Water 
Surfaces 

The term "set-up" may be defined to represent 
either the displacement of the water surface at the 
leeward end of the channel with respect to the undis- 
turbed level of the water or the difference of the 
displacements of the water levels at the windward and 
the leeward ends of the channel, that at the windward 
end being negative. The displacements referred to 
are the values that remain after correcting for the 
fall of the air pressure in the channel. In this 
investigation the term will be used in the second sense. 
Thus denoting the windward and leeward displace- 
ments by hi and h 2 , respectively, the set-up, S, is 



S=h2 — hi. 



(38) 



Observations of the set-up were made for depths of 
14.5, 11, 8, 6, and 4 cm with smooth water surfaces. 
In every observation of the set-up for a given wind 
velocity the following quantities were observed. 
First, the reading of the Venturi nozzle manometer, 
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which gave the average wind velocity in the channel, 
was noted. Second, the elevations of the water 
levels in the two extreme manometer wells were 
observed, the distance between the wells being the 
nearest possible to the fetch, that is, the length of the 
water surface acted on by the wind. The diameters 
of the manometer glass tubes being 2 cm, the meniscii 
were hardly affected by surface tension. The 
meniscus positions were read by a telescope against a 
centimeter scale with millimeter divisions, which was 
placed near and parallel to the tubes. The eleva- 
tions could be estimated with an accuracy of ±0.1 
mm. Third, the pressure fall in the channel air 
passage corresponding to the fetch was noted. This 
was done by connecting the extreme piezometer 
openings in the top of the channel to an inclined 
manometer, the pressure differences being read in 
heads of water. 

Because of the inclined solid bottom at the channel 
entrance the fetch varied with the depth of water. 
The manometer well distances, //, approximating 
the fetches were: 27 = 2,000 cm for 2Z = 14.5 cm; 
2/ = 1,850 cm for H=ll and 10 cm; and L'= 1,725 
cm for H=6 and 4 cm. The set-up values observed 
with these lengths were reduced, after correcting for 
the pressure fall in the air space, to a common fetch of 
Z,z=i ? 900 cm. As the curve of the surface displace- 
ments to a first approximation is linear, the reduction 
to a standard fetch is a matter of direct propor- 
tionality. That is, if S' is the set-up corresponding 




to L'j and if the difference L'—L is a small part of L, 
then 



Figure 7. Wind set-up in the absence of waves. 



S=S'L/L'. 



(39) 



The results of the set-up observations, determined 
in this manner, for various depths of water are shown 
in figure 7. The set-ups observed with smooth sur- 
faces will be denoted as S x . In the figure S u ex- 
pressed in centimeters, is plotted against V 2 , the 
square of the average wind velocity, where V is ex- 
pressed in centimeters per second. For a given 
depth of water and for a given fetch the set-up is 
proportional to the square of the wind velocity. 
From the data of figure 7 the values of Si for the 
different depths at a common wind velocity V= 1,000 
can be obtained. These are shown in figure 8. It 
is seen that the distribution of the data points is 
nearly linear, and the straight line drawn through 
the points has a slope of one to one. The interpre- 
tation, therefore, is that for a constant wind velocity 
and for a constant fetch the set-up is inversely pro- 
portional to the depth of water. Expressed mathe- 
matically, 

Si^V 2 /H. (40) 

Rendering the relation dimensionally correct we have, 
since the quantities g and L were constant in the 

tests, 



s 1 _ A v\ 



(41) 



where A is a numerical constant having the value 



yl=3.30X10- 6 . 



(42) 
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These relations give the law of wind set-up where 
the wind is not of sufficient strength to excite waves 
and therefore the water surface is smooth. 

3. Relation of Water Surface Velocities to Wind 
Velocities 

The argument regarding the development of the 
final formula of the preceding section is not com- 
plete. It remains to see whether the coefficient, A, 
of the wind velocity term is in every respect a nu- 
merical constant unaffected by the viscosity of the 
liquid. This aspect of the problem will be investi- 
gated later. For the present the inquiry will be 
facilitated by first establishing the relationship be- 
tween the surface velocities of the water and the 
wind velocities. 

An effort was made continually during the tests 
to measure the water surface velocity to see whether 
a relationship existed between this and the air 
velocity, as affected by the depth of liquid and the 
viscosity of liquid. It was also thought desirable 
to see if the relationship was modified by the fact that 
the liquid surface was smooth or covered by waves. 
The change in the viscosities was brought about 
invariabty by using sugar solutions of various con- 
centrations. 

To obtain the surface velocities, small paraffin 
particles were dropped from a perforation in the 
channel roof, and the movement of the particles 
along a given distance was timed by a stop watch. 
At the commencement of such observations, the 
velocities were measured over successive intervals 
covering the entire channel from the beginning to the 
end. It was found during numerous tests of this 
type that the surface velocities did not change with 
distance when the liquid depth was large and the 
viscosity of the liquid small. Only when the depths 
were small and the liquid viscosities large would the 
surface velocities vary with distance. Thus in the 
later experiments the surface velocity was measured 
over small lengths. 

The velocity values for a given wind velocity were 
based on the average value of three to five individual 
determinations. Usually the floating particles did 
not maintain a constant distance from the wall 
during their travel. The meander of the particles, 
however, was such that it could be taken to be 
equivalent to the motion of a hypothetical particle 



keeping a constant midway position between the 
channel axis and the glass wall. A study of the 
lateral distribution of the surface velocities did show, 
and the evidence will be reproduced later in this 
section, that the motions of such particles give very 
nearly the average value of the surface velocities 
taken in the lateral distribution. 

Representing the average surface velocity by 
u s , the question to be examined is the dependence of 
u s on the wind velocity F, on the liquid depth H f 
and on the kinematic viscosity v of the liquid. 
Obviously, the relationship is of the type 



V 



~K^> 



(43) 



For the establishment of the relationship the 
measurements obtained with water depths of 4, 6, 
8, 11, and 14.5 cm were first considered. In some 
of the tests the water was charged with soap or with 
detergent. These observations, however, did not 
yield Reynolds numbers of sufficiently low values. 
To supplement the data other observations were 
undertaken with liquid depths of 8, 4, and 2 cm, 
using heavy sugar solutions. 

The data accumulated from these extensive tests 
yielded the values of the parameters u s /V &nd u s H/i>, 
which are plotted in figure 9. Obviously, when the 
Reynolds number is small, the viscosity of the liquid 
has a very marked effect on the surface velocity. 
In the logarithmic plotting the distribution of the 
points of observation for low Reynolds numbers is 
linear, and the line drawn has the slope one in two. 
Accordingly the law for the surface velocities when 
the Reynolds number is small is 



v 
V 



=*(«£)», 



(44) 



where the coefficient K is a pure number having 
the value 



K=7.QX10~\ 



(45) 



Starting with a Reynolds number of approximately 
1,000 the ratio of the two velocities tends slowly to 
a constant value indicating that above this value 
the effect of kinematic viscosity is very slight. 
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Figure 9. Variation of water-surface 
velocities with Reynolds Number. 
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Judging from the curve drawn in the figure the 
limiting value of the ratio is 



V 



= 0.033. 



(46) 



In establishing the part of the graph above a 
Reynolds number of 1000, test data were used with- 
out regard to the condition of the liquid surface 
that happened to prevail. That is, some of the data 
are from runs where the surface was covered with 
waves and other data are from runs where the sur- 
face was smooth. Where waves were present, the 
surface velocities used were mean values from the 
crests and from the troughs. Due to the close 
grouping of points the inference is that the surface 
velocities are not affected by the presence of the 
waves. This surely was an unexpected result, and 
due to its importance the question was further 
investigated in two additional series of experiments. 
In one of the series the water of the channel was 
treated with detergent. During the measurement 
of the surface velocities the surfaces were smooth 
and calm. The data from this series of tests are 
shown in figure 10. For the second series the 
channel water was clean, and the observations were 
made with air velocities yielding substantial wave 
formations. The data from the latter series are 
shown in figure 11. In the figure are drawn two 
straight lines; the broken line represents the average 
values from the tests with water surface remaining 
calm; the continuous line represents the average 
values of the individual data shown in the figure. 
The disparity between the two curves is less than 
the scatter of the individual observations, and it 
may be stated on this basis that the presence of 
waves does not appreciably affect the relation 
between the surface velocities and the wind velocities. 
In the above discussion, the ratio of the two veloc- 
ities is treated as a function of the Reynolds number 
based on the surface velocities. It is equally feasible 
to take the Reynolds number based on wind veloci- 
ties. This is especially useful for comparison with 
prototype data. The data for this case are given in 
figure 12. Accordingly for Reynolds numbers below 
30,000 the ratio u s /V varies inversely with Reynolds 
number; for higher Reynolds numbers the ratio 
tends to a constant value of 0.033, roughly. 
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Figure 11. Water-surface velocities for large wind velocities 
in the presence of waves. 
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Figure 10. Water-surface velocities for large wind velocities 
in the absence of waves. 
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Before closing the matter of the present section, 
the evidence that the measured surface velocities 
equal the averages taken across the width of the 
channel must be shown. The establishment of the 
lateral distribution of the surface velocities proved 
to be a very difficult task. Much time was spent 
and the observations were repeated many times until 
some assurance about the results was obtained, 
working always with statistical averages. Without 
attempting a description of the methods used, the 
final results are shown in the form of graphs in 
figure 13. The lower group of points represents test 
data for the case where the liquid surface was smooth 
and free of waves. The data points are differently 
indicated and correspond to liquid depths varying 
from 2 to 8 cm and to velocities varying from 350 
to 700 cm/sec. As the points group so closely among 
themselves for a given distance from the walls, it 
did not seem necessary to indicate the velocities 
and the depths corresponding to specific symbols. 
The broken line gives the coordinate coinciding with 
the average velocity. According to these results a 
particle on the water surface at a distance of 2.2 cm 
from the nearest wall moves with a velocity equal 
to the average of the velocities in the lateral traverse. 
The upper group of points in the figure was obtained 
for a liquid surface covered with waves. The num- 
ber of observations was as large as for the points 
below; the observations being made with various 
depths and with various wind velocities. In order 
not to confuse the appearance of the graph, only 
averaged values are shown since no systematic varia- 
tions could be discerned in the data with respect to 
depth of liquid or to wind velocity. It appears as 
though the surface velocities vary less in a lateral 
traverse when waves are present. In this case also 



the surface velocity at a distance of about 2.2 cm 
away from the nearest wall represents practically 
the average velocity of a lateral traverse. 

4. General Law of the Surface Gradient in the 
Absence of Waves 

The regularity seen in the relationship of the sur- 
face velocities with the Reynolds numbers based 
either on the surface velocities or on the wind ve- 
locities suggests that perhaps a similar relationship 
should exist for the surface gradient. 

It was remarked in the theoretical analysis of 
purely viscous conditions that the surface gradient 
may be expressed as 



dh__6vu s 
d^~~gH~ 2 



(31a) 



This might be generalized to include the case of 
turbulent flow by writing 



where 



dh u 2 8 



* 



-*.(*?> 



(47) 

(48) 



The quantity 4> s may be referred to as the surface 
coefficient. Accordingly, the surface coefficient for 
the case of purely viscous flow is 



<i>s 



nil 



(49) 



The set-up observations that were made with water 
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Figure 14. Coefficient of surface 
slope in the absence of waves as a 
function of Reynolds Number. 



369 



charged with soap or detergent may be used for 
evaluating # s corresponding to large Reynolds num- 
bers, for in this case the surface gradient can be 
written also as 



dx L 



(50) 



The relation between the wind velocities and the 
surface velocities is given by the graph in figure 12. 
As regards <f) s for small Reynolds numbers, set-up 
observations were made with small depths (2 and 4 
cm) and with sugar solutions of various concentra- 
tions. In these latter tests, owing to the large 
viscosities of the solutions, waves could be excited 
only with winds of large velocities and therefore 
there was a range of velocities that left the liquid 
surface calm and smooth. As an illustration, the 
set-up observations made with a liquid of kinematic 
viscosity ^=0.111, a water depth H=2.0 cm, and 
a fetch L— 1,760 cm, are shown in figure 18. 

The determination of the surface gradient coefficient 
as a function of the Reynolds number is shown in 
figure 14. It is apparent for low Reynolds numbers, 
that if the coefficient be defined as in eq 47, its value 
is affected very markedly by the Reynolds number. 
The distribution of the data points is sufficiently 
linear to allow the drawing of a straight line through 
them. The line actually drawn conforms to the 
theoretical expression given by eq 49, and since the 
line represents the positions of the observed points 
quite accurately the theory is verified satisfactorily. 
There was another deduction of the theory that 
appeared from limited observations to be verified. 
According to eq 29 the particle velocity is zero at a 
point two-thirds of the total depth measured from 
the bottom. Observations of small particles that 
had accumulated in the channel for one reason 
or another indicated this deduction to be sub- 
stantially true. 

For high Reynolds numbers, as may be observed 
from the data of figure 14, the effect of Reynolds 
number is small and the coefficient tends to a con- 
stant value. The deduction from the data presented 
is that the surface velocities define completely the 
regime of flow in the liquid below, the flow assuming 
either the character of a viscous flow or a turbulent 
flow. The passage from one type of flow to the 
other type does not show the sharp and sudden 
transition ordinarily seen in fluid motion phenomena. 

For practical applications it is desired to define 
the coefficient of the surface gradient in terms of the 
wind velocities instead of the liquid surface velocities. 
In this new representation 



dh 



V 2 



and 



dx~~* a gH 7 



<t>a = <l>a 



(?) 



(51) 
(52) 



where <j> a is the coefficient of the surface gradient in 
relation to wind velocity. The transfer from <£ s to <j> a 



is readily made by using the data in figures 9 and 14. 
The determinations of €/> a thus made are shown in 
figure 15. The surface gradient coefficient <j> a is 
independent of the Reynolds number VH/v. 
Specifically 

g=3.35X10- If < 53 > 

This value of </> a =3.35X10" 6 is to be compared 
with the value of ^L=3.30X10~ 6 appearing in eq42. 
This then completes the argument that the coeffi- 
cient, A, appearing in eq 41 is actually a numerical 
constant and in particular is independent of liquid 
viscosity. 

The result shown above is one of remarkable 
simplicity. In terms of the shear acting on the 
water surface, the gradient is 



dh nr s 

dx pgH' 



(22a) 



and comparison with eq 53, using the average value 
of <j> a from eq 42 and eq 53 yields, 



7ir s =:3.32X10- 6 p^ 2 



(54) 



where p is the density of the liquid. Replacing p 
with the moist air density p a by using the relation, 



Pa =1.169X10- 3 p, 
the expression for the surface stress is 
nr s = 0.002Sp a V 2 . 



(55) 



(56) 



If the value 1.5 be assigned to n, the value of the 
surface stress, using the ordinary form of the formula 
for such stresses, is 



= 0.0037 



(57) 



If it is supposed that the height of the air passage 
in the present experiments has no effect on the stresses 
developed on the surface of the liquid, then the shear 
obtained might be compared with the shear evalu- 
ated according to some boundary layer theory ap- 
plicable to the case of a moving liquid surface. 
Such a theory is not yet available, and the next best 
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Figure 15. Coefficient of surface slope in the absence of waves 
as a function of Reynolds Number based on wind velocities. 
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choice is therefore to compare our results with the 
well-known theory of the resistance of smooth plates. 
In terms of the drag coefficient C d , 



where 






(58) 
(59) 



and v a is the kinematic viscosity of air. In the 
tests, the results of which are now being examined, 
the Reynolds number VL/v a ranged from 10 6 to 
10 7 . The drag coefficient for this range as deter- 
mined by Gebers [5] and by Kempf and Kloess [6] is 



<7 d =0.0032, 



(60) 



and this agrees very well with the corresponding 
coefficient in eq 57. 

It would thus appear that the stress developed 
over the smooth surface of water under the action of 
flowing air could be determined with sufficient 
accuracy by using the ordinary theory of smooth 
plate resistance. 

For emphasis, the foregoing results will be sum- 
marized. Shear acting on the liquid surface in the 
absence of waves is independent of the depth of 
liquid and the kinematic viscosity of the liquid. The 
magnitude of the shear is almost identical with the 
shear that may be associated with the skin friction 
of smooth plates. In this comparison one may 
assume that the ratios of the surface shear to the 
bottom shear is two; that is, n=1.5. The value of 
the ratio is exact in the case of viscous flow; it is not 
known precisely for the case of turbulent flow. The 
prevailing flow in the liquid below the surface is 
determined solely on the basis of the surface veloci- 
ties, the depth, and the kinematic viscosity of the 
liquid. For purely viscous flow, the flow is in 
accordance with the simple one-dimensional theory. 
The exact conditions and the details of the mecha- 
nism of flow in the case of turbulent flow are not 
known. These could be made the object of a fruitful 
and a promising special study. 

V. Wind Set-Up in the Presence of Waves 
1. Effective Air Velocities 

In the observations of set-up with smooth water 
surfaces the variations in the value of the air velocities 
with distance along the channel axis proved to be 
small, and they were ignored. On the other hand 
the corresponding variations were large when the 
surface of the liquid was covered with waves and due 
account was taken during the reduction of data. 
Two causes were operative in producing the varia- 
tions of the air velocities with distance. First, waves 
induced comparatively large stresses over the water 
surface and thus a greater rise in water level in the 
leeward end of the channel. Second, as the heights 



of the waves increased with distance, the elevated 
crests caused a further decrease in the height of the 
air passage. The air velocity variations with 
distance were especially marked for the larger depths 
of water because of the small depth of the channel. 

Under these conditions the air velocity values 
derived from the Venturi nozzle are only nominal 
velocities, and these had to be modified or corrected 
to obtain the effective air velocities. The basis of the 
corrections is as follows: 

The local surface gradient at the point x may be 
written as 



dh n V 2 
ax (jH 



(61) 



where the coefficient is 



c=c (jL , w , L'~r)' (62) 

which is an assumption of the most general applica- 
tion. If V be constant, the integration of eq 61 
yields 

V 2 C L V 2 

S=h 2 -h 1 =^ T f\ Cdx^C m L^ (63) 
gHJo gti 

Supposing that V is not constant but V=V n -\-AV y 
where AV is a small fraction of V n) the same integra- 
tion gives 

or expressed in a different manner 

S=^] [ C m L+2r ] C m £ X ^dx], (64) 

where rj is a quantity larger than unity. Comparing 
with eq 63, 

which gives the effective wind velocity, V, in terms 
of the nominal wind velocity, V n > What the exact 
value of the coefficient rj should be is questionable. 
We have assumed 77 = 1.5 as a plausible value for 
present applications. 

To apply the correction to those measurements of 
set-up that showed considerable restrictions in the 
air space above the waves it became necessary to 
observe the variations of AV/V n along the channel 
axis for the different depths of water and for the 
different air velocities. To give an idea of the 
variations noted, the data for a depth of 11 cm are 
shown in figure 16. Using these data and also 
similar data for other depths, the relationship be- 
tween the effective wind velocity and the nominal 
velocity was established as shown in figure 17. 
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axis. Water depth= 1 1 .4 cm. 



0.12 



SYMBOL 


H 




cm 


O 


14.5 


□ 


11.0 


V 


8.0 


o 


6.0 


A 


4.0 



c .06 




500 600 700 600 900 1000 

V n , Cm /sec. 

Figure 17. Effective wind velocities in terms of the nominal 
wind velocities. 



2.0 



r 



H = 2.0 cm 

l/= O.I 110 cm /S«C 

L = 1760 cm 





60 



80 



, Cm %* 



Figure 18. The partitioning of set-up when waves are present. 
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2. Dependence of Set-Up on Depth of Water When 
Waves Are Present 

Great difficulty was experienced in the beginning 
when an attempt was made to establish the relation- 
ship between set-up and the depth of water for the 
case where waves were present. The difficulty was 
overcome partly by supposing that it is permissible 
to break the observed set-up into two parts, Si and 
S 2 , the first being due to skin friction alone and the 
second part due to form resistance of the wave crests. 
It was supposed further that the magnitude of S x 
is not affected by the presence of waves, and indeed 
that it would be the value that would have been 
obtained for the prevailing air velocities with a smooth 
water surface. Thus, if S is the observed total 
set-up, the part S 2 due to the waves is simply 



where 



S 2 =S — Si, 
Si = 3.30XlO~ Q ~L. 



(66) 
(41a) 



The latter relationship was fully discussed in the 
previous sections. 

At the moment the dynamical significance of the 
procedure of partitioning is not clear. The only 
justification that may be mentioned now is the 
consistency and the reasonableness of the results 
obtained. The idea was suggested by the fact 
that the relationship between the water surface 
velocities and the wind velocities is not affected by 
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Figure 19. Partial set-up due to waves in tests with small 
water depths. 



the presence of waves. Eventually, of course, we 
wish to consider the dynamical aspects of the 
phenomenon of total resistance in a future investiga- 
tion. For the present, however, our main interest 
is in the parametric representation of measured 
quantities and the validity of these representations 
for the quantitative evaluations of the set-up. 

As an example of partitioning the set-up into parts, 
the data in figure 18 are presented. The data were 
obtained with a sugar solution of kinematic viscosity 
y=0. Ill cm 2 /sec. The surface length was Z= 1760 
cm, and the depth of water was H=2 cm. Due to 
the high value of viscosity, the surface remained 
smooth for an appreciably large range of wind veloci- 
ties, as can be seen from the curve. For this range 
the set-up varies with the square of the wind velocity 
as it should. Waves were initiated when the wind 
velocity exceeded 650 cm/sec. The curve indicates 
that the increase in the value of the set-up starts 
with a jump. The projection of the straight line 
passing through the points of observation correspond- 
ing to small wind velocities is indicated by dashes. 
The quantity S 2 shown in the diagram, which is the 
excess over the dashed line, is taken to represent that 
part of the set-up due to waves. 

The same method of graphic partitioning was em- 
ployed in the results of the tests made with water 
depths of 14.5, 11, 8, 6, and 4 cm. As the tap water 
used in the tests was relatively pure, waves of sub- 
stantial size were obtained. The data are given in 
figures 19 and 20. The straight lines represent the 
magnitudes of S u whereas the curves drawn through 
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the symbols show 'the excess of the observed set-up 
with respect to Si. These two quantities, Si and S 2 , 
are expressed as functions of the squares of the wind 
velocities. As is seen from the figures, S 2 commences 
to increase after V has exceeded a certain critical 
wind velocity, V . The figures suggest that the in- 
crease starts with a jump. This fact will be ignored, 
and it will be assumed that the increase is from a null 
value and the changes are gradual. The curves rep- 
resenting S 2 are drawn accordingly. 

On the basis of dimensional reasoning it is to be 
expected that 



& =/ f_V_ (£^ 3 P_ H\ 
L Jl \J0 y'gH'pa'Lf 



and as to the identification of V , 

■yfgH j2 \^gH'pa'Lj 



(67) 



(68) 



Eliminating the parameter containing the kinematic 
viscosity, 



L 



3 \jgFI ^/gH Pa LJ 



(69) 



Now, if we take V as 350, 360, 380, 420, and 440 
cm/sec, corresponding to water depths of 14.5, 11,8, 
6, and 4 cm, respectively, then the additional set-up 
S 2 due to the waves may be represented as 



s 2 =c(v- Vo) 2 



(70) 



The evidence for the statement is shown in figure 21. 
Since in this group of data the quantities i, g, /x are 
not changed, the factor C may involve the depth H. 
To determine the nature of dependence on the depth, 
the values of S 2 corresponding to (T 7 — V ) 2 =50X 10 4 
are first read from the lines in figure 2 1 and are then 
transferred to the logarithmic plot of figure 22. The 
disposition of the points definitely lacks the even 
regularity of alinement that one naturally wishes to 
find. The points may be approximated by a straight 
line ; but there is considerable latitude in the selection 
of the slope of the line. If one adopts the straight 
line drawn in the figure, then the equation for the 
partial set-up due to waves is 



S* 



=B 



(V-Vo) 2 



fH\' 2 



(71) 



where 7?=2.08X10" 4 , and T r has the values men- 
tioned above. If we disregard the small variations 
of V Q for the different depths, the mean T 7 o=390 
cm/sec is a suitable value. Adding to this the part 
due to skin friction, the expression for the set-up in 
the presence of waves is 



8 



gH gH 



(?)"• 



(72) 



where 



A=33.0X10" 6 , 5=2,08X10" 



It is understandable that the expression represents 
with sufficient accuracy the data that were observed. 
Yet, neither the reality nor the dynamical significance 
of the quantity V appearing in the formula can be 
discussed at the moment in a manner satisfactory to 
us. For this reason the quantity T /T will be referred 
to as the formula characteristic velocity. 

3. Effect of Viscosity on Partial Set-Up Due To Waves 

It was shown previously, working with the data 
for smooth water surfaces, that the partial set-up 
Si is independent of viscosity. It now remains to 
investigate what effect, if any, viscosity will have on 
the partial set-up due to waves, that is on S 2 . 

To consider this matter, a series of tests was car- 
ried out with a constant liquid depth of 8 cm. In 
addition to water, sugar solutions of four different 
concentrations were used, thus making altogether 
five experimental liquids of unequal kinematic vis- 
cosities. The setups observed with these liquids 
for an effective surface length of 1,900 cm are shown 
in figure 23. The two components, Si and S 2 , of the 
total setup are plotted as functions of the squares 
of the wind velocities. It is to be noted that the 
Si line is common to all the liquids and is independent 
of viscosity. After drawing the curves through the 
S 2 points of observation, values of S 2 were read from 
these curves, the square roots of which are plotted 
in figure 24 against the wind velocity. The disposi- 
tion of the points is linear. This would mean that 
the form of the expression in eq 70 is valid for 
liquids of varying viscosity. The lines have practi- 
cally the same inclination independent of viscosity, 
and thus the value of C in eq 70, and hence the value 
of B appearing in eq 71, is not affected by viscosity. 
The intersections of the lines with the axis of ve- 
locities give the formula characteristic velocity, V , 
which varies with the viscosity. Specifically, cor- 
responding to the successive viscosities y= 0.0092, 
0.0154, 0.0217, 0.0351, 0.0600 cm 2 /sec, the formula 
velocities are T r =380, 630, 680, 720, 900 cm/sec, 
respectively. 

Summarizing, the only effect of the increased 
viscosities is to increase the formula characteristic 
velocity. As a result, then, the total setup observed 
for a given wind velocity decreases when the kine- 
matic viscosity of the liquid is increased. 

4. An Analysis of Hellstrom's Data on Setup 

In deriving the law of setup expressed by eq 72, 
the assumption was made tacitly that the limited 
height of the air passage over the surface of the 
water in these experiments does not substantially 
modify the magnitude of the resulting setups. 
Had it not been for the fixed construction of the 
channel, this circumstance would have been studied 
experimentally by varying the height of the air 
passage for a given depth of water. Fortunately 
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Figure 21. The relation of the formula characteristic velocity 
to the partial set-up due to waves. 
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Figure 23. Effect of viscosity on partial set-up due to waves- 



SYMBOL V 




c ™%ec 


U 


.0092 


O 


.0154 


A 


.0217 


o 


.0351 


V 


.0600 



l£ 




800 1000 

V cm / 
v * /sec 



1400 



Figure 24. The determination of formula characteristic 
velocities. 



375 



the experiments of Hellstrom are of value in discuss- 
ing the effect of size of air passage by comparison 
of results, as the channel used by Hellstrom had an 
air passage of relatively large depth. 

The data in figure 25 give the relationship between 
surface shear and depth of water for various wind 
velocities as found by Hellstrom. One modifica- 
tion was made on our part; that is, to conform to 
our results channel central velocities were changed 
to mean velocities in the channel cross section. The 
data in this figure may be made the basis on which 
to deduce the set-up magnitudes in the Hellstrom 
channel. The channel had a length of 300 cm and a 
depth of 50 cm. The average depth of the water 
in the experiments appeared to have been varied 
between the limits of approximately 2 and 4 cm. 

From the relation 



dh_ nr s 
dx pgH 



(22b) 



if the set-up is small, can be derived the expression 

JlToL 



S = 



pgB*. 



(73) 



which connects the set-up with the effective average 
wind force over the entire fetch. Since Hellstrom 
has assigned to n the value 1.50, the proper relation 
to be used in deducing the set-up in his experiments is 



S = 



1.5 r s L 



>gH 



(74) 



Taking the data from figure 25, sequences of set-up 
values were obtained for the depths of 2, 3, and 4 
cm. These computed values are given by the points 
plotted in figure 26. 

To compare these with the results of the present 
investigation it suffices to determine similar values 
of set-up, using eq 72. Now, it will be remembered 
that for water the formula characteristic velocity 
showed slight variations from the mean, T7 =390 
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Figuke 25. The relation between wind velocity and depth of 
water in Hellstrom's experiments. 



cm/sec, as the water depths where changed. We wish 
to assign a similar tendency to be valid in the Hell- 
strom tests. Accordingly we assign to V the mag- 
nitudes 410, 390, and 370 cm/sec to apply to the 
depths of 2, 3, and 4 cm, respectively. The com- 
puted values of the set-ups on this basis, using our 
general formula, are shown likewise in figure 26. 
The curves of the computed set-up values form a 
band, and the Hellstrom values are within this band. 
Moreover, the sequence of the curves is the same as 
the sequence of the Hellstrom set-up values as the 
depths are changed. This shows that there is quite 
a satisfactory agreement between the general formula 
and the results from the Hellstrom channel. Ac- 
cordingly, we are permitted to infer that the restric- 
tive effect of the small air passage in our experiments, 
if present, must have been small indeed. 



VII. The Probable Meaning of the Formula 
Characteristic Velocity 

1. Dimensional Requirements of the Monomial Set-Up 
Formulas 

Before we consider various topics related to the 
question of the formula characteristic velocities, it 
becomes necessary to see if the total set-up in asso- 
ciation with waves can admit representation by a 
single term. The general formula developed in this 
investigation is at variance with this possibility, 
whereas a monomial formula is implied in the work 
of Hellstrom. 

Hellstrom proposes that the relation between the 
wind force and the wind velocity be of the form 



-HV 3 



(75) 



In view of the relation 
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Figure 26. Comparison of set-ups from the Hellstrom experi- 
ments with values obtained from equation 72. 
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dh 7iT s 

dx pgH' 



(22c) 



one must now infer that the corresponding expression 
for set-up is 

j^A'V*, (76) 

where A' is independent of the depth of water. 

The set-up values that were deduced previously 
and considered during the analysis of the Hellstrom 
data are reproduced in figure 27. There are repro- 
duced also the set-up values observed in the tests of 
the present investigation with water for the region 
of high wind velocities. The graph representation 
is logarithmic and the observation points of the two 
groups of tests, one from Hcllstrom's channel which 
was 300 cm long and the other from our channel, 
which was 1,900 cm long, aline themselves linearly. 
In both cases the inclination of the lines representing 
the points may be taken to be approximately 3 to 1. 
Thus the formula of eq 76, applies equally well to 
the two groups of data from water surfaces covered 
with waves. 

Now in order that eq 76 may represent the physical 
law of set-up in laboratory channels, the right-hand 
member must be a dimensionless expression since 
the left-hand member is a dimensionless ratio. It will 
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Figure 27. Monomial formula of set-up with waves. 
i>= 0.009 to 0.010 cmVsec. 
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be remembered that the data now being considered 
are from two groups of tests where the kinematic 
viscosities were the same. Accordingly, the general 
law could be 



§=-(yfzK 



fgLI> 



(77) 



The right member contains a factor involving the 
kinematic viscosity. The influence of this factor 
may be determined by considering the set-up data 
obtained with highly viscous sugar solutions subjected 
to large wind velocities. These data are shown in 
a separate logarithmic plotting in figure 28. The 
distribution of the points for a given viscosity ^ are 
linear, and furthermore the inclination of the lines 
are in agreement with the expression in eq 77. 

If the S values for the velocities V= 1,400, 1,000, 
and 800 cm/sec be considered, as in figure 29, it 
is seen that these may be approximated by straight 
lines with a slope of — }{. Accordingly, if the expres- 
sion in eq 76 does represent a physical law, it must 
be of the form 






) 3 (^T 



(78) 



It now remains to be seen if this is actually the 
physical law for the set-ups observed in the laboratory 
channels. Owing to the manner of the derivation, 
the separate exponents represent the effects of veloci- 
ty and kinematic viscosity. To examine the con- 
tingency that these exponents represent likewise 
the effect of channel length, the last expression will 
be written in the form 



S=A' 



V" 



04/ 3j/ l/3 



(79) 



Accordingly, wind velocities and kinematic viscosi- 
ties remaining the same, the set-up observed is 
independent of the length of channel. As this deduc- 
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Figure 30. Variation of initial fetch for the start of waves with 
wind velocity. 



tion is contrary to observation, it is to be inferred that 
the monomial form for expressing the set-up is 
invalid when set-ups are observed in laboratory chan- 
nels with the wind velocities sufficiently high to 
produce waves. 

2. Critical Wind Velocity for Wave Formation 

It was noted during the preparation of data for 
analysis that almost invariably the formula charac- 
teristic velocity exceeded the critical wind velocity 
for generating waves. It may prove to be of some 
significance if a relationship can be established be- 
tween these two velocities. This in turn requires 
that the law giving the critical wind velocity for 
wave formation be first established. 

At very low wind velocities the surface of the 
water remained calm. In fact if one studied the 
images of the vertical and the horizontal bars of 
windows reflected from the water surface, they 
appeared as even bands extending laterally and 
longitudinally over the water surface. On increas- 
ing the wind velocity, a tremor of the water surface 
was noticed, and the reflections of the bars now had 
the appearance of sinusoidal bands. This condition 
extended from one end of the channel to the other, 
except that there was a tendency for the amplitudes 
of the irregularities to increase somewhat with 
distance. With a further increase of the wind ve- 
locity the tremor of the water surface was consider- 
ably augmented until finally the reflected shadows 
were confused. Experience showed that this state 
of the water surface was precedent to the appearance 
of waves, and in fact only a slight further increase 
in the wind velocity would be sufficient for the forma- 
tion of waves. 

The waves first formed were two dimensional, 
their length being estimated to be close to 4 cm. 
The wave heights were less than 2 mm, but this 
value could not be definitely determined since it was 
effected by the meniscus at the glass wall. Initially 
when the waves were formed, they did not cover the 
entire area of the water surface but were seen first in 
the leeward part of the channel. This would mean 
that an initial fetch of length, X w , is necessary for 
the waves to appear, or more precisely to be visible 
to the eye. With increasing wind velocities the 
length of the initial fetch gradually decreased. To 
examine the magnitude of the decrease a few observa- 
tions were made at water depths of 14.5 and 8 cm. 
The data from these observations are shown in figure 
30. It is obvious that for these observations the 
initial fetch is inversely proportional to the square 
root of the velocity of the wind. When the initial 
fetch decreased to a length of about 1 or 1% m, a 
new condition of the surface was manifested. The 
area between the windward water edge and the 
locality where the initial two-dimensional waves of 
the wavelength of 4 cm were first noted was covered 
with ripples. These ripples were three-dimensional 
and were about 2 cm long. One noticed with 
amazement the sharp and distinct issuing of the 
initial two-dimensional waves of wavelength of 4 
cm at the lower end of the area of the ripples. With 
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greater wind velocities the initial area of the ripples 
was somewhat decreased; but in these investigations 
the length of the area never fell below 60 cm. 

When the initial fetch for the startof the waves was 
small, which occurred at high wind velocities, wave 
height and wavelength increased with distance. 
With some exceptions, the increase of the wave 
height was proportional to the square root of distance, 
whereas the increase of wavelength was proportional 
to the distance. When the initial fetch was large, 
which corresponded to small velocities, it was ob- 
served at times that the waves once formed would be 
damped out, but reappeared near the leeward end of 
the channel. 

The condition of variable fetch at the start of the 
waves and the tendency of waves to damp out when 
the wind velocity was small made it difficult to de- 
termine the particular wind velocity to be associated 
with the formation of waves. The matter of obser- 
vation of the critical wind velocities for wave genesis 
was resolved as follows. Three initial fetches, of 
lengths 1,000, 800, and 600 cm were selected and 
the corresponding wind velocities determined. The 
average value of the three velocities was taken to 
represent the critical wind velocity, V c , for the wave 
formation. When using water alone, the critical 
wind velocities were determined for depths of 14.5, 
11, 8, 6, and 4 cm. When using more viscous solu- 
tions, the depths were 8, 4, and 2 cm. 

In the presentation of data we follow the procedure 
of dimensional analysis. If it be assumed that the 
critical wind velocity is independent of fetch, then 



<*v 














































/ 


/ 






















/ 


/ 

o 


o 
















o 




/ 


/ 

o 


' 


o 






























Q 





































s° 
















[o» 6 








3 / 


/ ° 




















o 




















>° 






o/ 


1 
















^ 


V c =I6.2Qji/^] 3 






























































1 



























O.I 



0.? , 0.4 



0.6 0.8 1.0 



Figuke 31. Critical wind velocity for the start of waves. 



the dependence can be on the depth of water, the 
kinematic viscosity of the liquid, and the densities 
of the air and the liquid. Accordingly, 



V e 



gH 



= J([9VP/Pa\ 



(80) 



Conforming with this, the observed data are pre- 
sented in the logarithmic plotting in figure 31. De- 
spite the difficulty encountered in the observation 
of the critical wind velocities, the distribution of the 
points is sufficiently regular to allow the drawing of 
a straight line to represent them. Since the slope 
of the line drawn is unity, the critical wind velocity 
is independent of depth of liquid. Accordingly the 
desired form of the law is, very simply, 



V c =l§.2(g P vl Pa yi\ 



(81) 



Interestingly enough the form of the criterion 
accords well with the form that Harold Jeffreys [7] 
has derived theoretically. The only disagreement 
is in the numerical values of observation. The 
observations made by Jeffreys in open lakes and 
rivers suggest that the numerical value of the coeffi- 
cient is about one-third as great as that given in eq 
81. The cause of the disagreement may be the 
possibility that the critical value in reality depends 
also on the fetch of the initial waves noted. This 
aspect of the question deserves an examination in 
channels of variable lengths. 

3. Correlation Between the Critical Wind Velocity 
and the Formula Characteristic velocity 

As mentioned previously, it was noted during the 
work of data reduction that a correlation may exist 
between the critical wind velocity and the formula 
characteristic velocity. Since, as shown above, the 
factors affecting the critical wind velocity assume a 
very simple relationship, it would likewise seem 
sufficient to consider the dependence of the formula 
characteristic velocity on the kinematic viscosity 
and the channel length in the manner carried out in 
figure 32. The open circles in the figure are the 
values obtained during the analysis of the set-up 
due to waves, and the values of the associated 
formula characteristic velocity were indicated pre- 
viously. The values corresponding to the black 
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Figure 32. Variation of the formula characteristic velocity with 
viscosity. 
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circles were obtained in a different manner. The 
method is suggested from the graph in figure 18. 
The velocity concurrent with the sudden increase 
in the set-up values due to the presence of waves is 
interpreted as a value close to V . The full circles 
are from tests with small liquid depths and with high 
viscosities. 

The distribution of the points suggests that it 
should suffice to represent the data by a horizontal 
line. Thus the formula characteristic velocity is 
independent of channel length and is proportional 
to critical wind velocity for the generation of waves. 
More specifically 



v =2i.o(gpv/p a ) 



1/3 



(82) 



The reader unquestionably will suspect, and this 
very justly, that we have tended to an over-simplifi- 
cation in this representation. Although the use of 
the formula characteristic velocity has caused a 
simpler representation of the set-up relationship 
in the general formula, this velocity by itself must 
have reference to various hydrodynamic elements 
involving dissipation and relative velocities such as 
the velocity of the waves with respect to the velocity 
of wind and also the velocity of the drift current at 
the surface of the water. At the moment the 
author is not prepared to discuss these matters 
satisfactorily. If dissipation really plays a large 
part in forming the formula characteristic velocity, 
then the opinion may be expressed that with greater 
liquid depths the magnitude of the formula char- 
acteristic velocity decreases. 

4. A Possible Mathematical Basis for the Law of 
Partial Set-up Due To Waves 

As mentioned before, between the points rep- 
resenting S 2 values and the curves intended to 
represent the course of the data, there exist small 
disparities. This fact no doubt has a bearing on the 
accuracy of the determinations, first of the formula 
velocity, V , and second of the exponent of H/L 
occurring in the formula for the partial set-up S 2 , that 
is in the expression of eq 71. Accordingly, if one can 
formulate a mathematical basis that will justify the 
form of the law for S 2f independent of the actual 
methods of reductions used, this should prove to be 
very illuminating. 

Pursuant to our original ideas and opinions we will 
suppose that when waves are present each crest 
contributes an additional resistance, and the resist- 
ance, in its major part at least, is due to form 
resistance. Let R be the resistance. The length of 
waves being X, the added shear due to waves is 

r s2 =R/\. (83) 

As for the law relating to R we may suppose that 

R=cp a (V-Vo) 2 a, (84) 

where a is the half value of wave height. At the 



moment nothing will be said as to the cause or the 
origin of V Q . Thus, the added shear is 



T s2 =cp a (V— V ) 2 a[\. 
Substituting in the expression 



'dh\ 



Ts2 



(85) 



(22c) 



and integrating with respect to x, there is obtained 



S 2 



(V- Vo ) 2 { L a 



dx. 



(86) 



If the variation of aj\ with respect to x be known for 
the waves generated in shallow channels by the 
action of wind, the integration required above will 
lead to an expression for the set-up #2 due to waves 
only. 

Next it shall be supposed that the law of growth of 
waves in shallow waters is virtually 



and 



a^V^Hx 1 ' 2 , (87) 

\~V w H"*x> (88) 

for sufficiently large values of V. Hence, we have 



-(f) 



(89) 



That is, the steepness of wind waves in shallow water 
decreases with distance in the leeward direction. 
Since 



£ 



■dz^CzH 1 ' 2 !* 1 ' 2 , 



there results, by substituting in eq 86 
L gH \LJ 



(90) 



(91) 



which in fact is the form that was adopted in the 
analysis of the set-up due to waves. 

In the above derivation three separate assumptions 
are used. These are represented by eq 85, 87, and 
88. Each of these assumptions is amenable to ex- 
perimental examination. This work is now under 
way. 

VIII. An Application to Natural Conditions 

In his fine treatise Hellstrom gives some data on 
the wind effects of storms on Lake Erie. It will 
be interesting to compare the magnitudes of the 
set-ups occurring during these storms with the com- 
puted values to be obtained with the general formula 
of the present investigation. The observed set-ups 
at Lake Erie are shown in figure 33 by the plotted 
points. 
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Figure 33. Comparison of observed set-up in Lake Erie during 
storms with equation 87. 



In applying the general formula, eq 72, a statement 
is necessary regarding the value of the formula char- 
acteristic velocity appropriate for natural conditions. 
As was mentioned previously, the critical wind ve- 
locity for the genesis of waves on large bodies of 
water is about one-third as great as the correspond- 
ing value obtained in our laboratory channel. 
Furthermore, it was also shown that the formula 
characteristic velocity is about 1.3 times larger than 
the critical wind velocity. Thus it is to be expected, 
as the formula characteristic velocity is of the order 
of 100 cm/sec for large bodies of water, that only a 
minor error will be made if this velocity is omitted. 
With this understanding the set-up expression for 
natural conditions may be given in the form 



S 



= 3.30X10" 6 1 + 63 



.[, 






(87) 



and this should apply to large bodies of water, the 
confinement approximating the shape of a rectangular 
channel of a uniform cross section. 

The average depth of Lake Erie being 17.7 m and 
the length 377 km, computations on the basis of the 
above formula gave a set of values of the expected 
set-up shown by the heavy straight line in figure 33. 
It will be remembered that the difference between 
the surface elevations at the extreme ends of the 
lake is defined as the set-up. 

The agreement between the laboratory results and 
the observed set-up during storms may be deemed 
satisfactory. However, this relationship can neither 
be elaborated on nor discussed with a sense of fi- 
nality. In order to make a definite evaluation of 



the rigor of the general formula proposed on the 
basis of the tests of the present investigation, it is 
necessary to have reliable observational data from 
numerous lakes or reservoirs with different depth/ 
length ratios. Such data are not known to us at 
the present. 

IX. The Need of New Experiments 

One of the most significant deductions from the 
general set-up formula developed from these tests is 
in relation to the effective stresses at the surface of 
the water. Ordinarily it is supposed that the stress 
or the wind force is independent of the length of 
channel. The tests of the present investigation, 
however, show that the stress depends on the length 
of the channel. If this result holds in general, model 
experimentation on wind effects to simulate proto- 
type conditions and effects will be out of the 
question. 

Regarding this effect of channel length, two pos- 
sibilities may be imagined. First, the effect of the 
length on the effective magnitude of the wind force 
may continue indefinitely when the channel length 
is increased keeping the water depth constant. 
Second, the clfect may become gradually less impor- 
tant with increased length and less than that required 
by the general formula now given. 

These matters can be studied within certain 
limitations in a new series of tests with channels much 
longer than the channel used in the present investiga- 
tion. At the National Hydraulic Laboratory there 
are now available channels of sufficient length, and 
it is our intention to resume experiments when oppor- 
tunity permits. 



The author expresses his appreciation for the 
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